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Abstract. The Fermi Gamma-ray Space Telescope (formerly known as Gamma-ray Large 
Area Space Telescope, GLAST) was successfully launched on June 11 2008. Its main instrument 
is the Large Area Telescope (LAT), which detects gamma rays from 20 MeV to more than 300 
GeV. It is a pair-conversion telescope with 16 identical towers (tracker and calorimeter), covered 
by an anti-coincidence detector to reject charged particles. The calorimeter is a hodoscopic 
array of CsI(Tl) crystals, arranged in 8 alternating orthogonal layers, with a total thickness of 
8.6 radiation lengths. In this paper we will present the performance of the LAT, with special 
O ' attention to the calorimeter, which provides a good energy measurement up to 3 TeV. We will 

also review some of its scientific results after 4 years of operation, focusing on measurements 
which extend up to very high energy, such as the spectrum of the extragalactic diffuse emission, 



o 



the spectrum of cosmic electrons and the positron fraction. 
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1. Introduction 

' The main scientific goal of the Fermi Large Area Telescope [JJ is to study the gamma-ray 

sources from 20 MeV up to more than 300 GeV. There is a large variety of such gamma-ray 
sources, from galactic sources (pulsars, pulsar wind nebulae, supernova remnants, microquasars) 
to extragalactic sources (active galactic nuclei and gamma-ray bursts). After 3 years of operation, 
the LAT had detected more than 1800 gamma-ray sources. On top of individual sources, we 
also seek to measure the galactic diffuse emission and the extragalactic diffuse emission. 

The point source sensitivity is mainly the result of two characteristics of the telescope: its 
angular resolution and its total acceptance (the effective area times the field of view). The 
angular resolution drives the ability of finding a source on top of a relatively flat background. 
The total acceptance determines the number of gamma rays that can be detected from a given 
source. A large effective area is important, but a large field of view is also very important: it 
allows the telescope to observe a large portion of the sky at the same time, so that the total 
amount of time during which a given source is seen, is a large fraction of the total observational 
time of the telescope. Achieving the best angular resolution and the largest total acceptance 
were the main goals of the design of the LAT. 

Regarding the energy measurement, the LAT was optimized to have the largest energy range: 
down to about 20 MeV and up to 300 GeV and more. In this energy range, gamma-ray 
sources generally have a smooth power-law energy spectrum with a photon index about -2. 
Thus measuring their energy spectra only requires a moderate energy resolution. But, since the 
spectra fall strongly with energy, any significant migration of reconstructed energy toward higher 
energies would easily distort the energy spectrum measurement. So the energy reconstruction 



must minimize the risks of overestimating the energy. Though most of the gamma-ray sources 
have a smooth energy spectrum, some, like pulsars and active galactic nuclei, have energy 
cutoffs or breaks that we want to measure as precisely as possible. We also want to reach 
the best sensitivity to gamma-ray lines that could result from dark matter annihilation. As a 
consequence, achieving the best energy resolution is also important for the LAT. 

2. From EGRET to the LAT 

In order to detect gamma rays above 10 MeV and measure their directions and energies, a space- 
based detector needs three sub-detectors: a tracker, in which gamma rays convert to pairs, to 
measure their directions, a calorimeter to measure their energies and an anti-coincidence detector 
(ACD) to reject charged particles. The tracker is made of some passive material in which the 
conversion occurs and some active material that tracks the e + e~ pair. The effective area depends 
on the fraction of converted gamma rays, which increases with the amount of converter. But 
below ~ 10 GeV, the angular resolution is limited by multiple scattering, which also increases 
with the amount of converter. So a compromise must be found between increasing the effective 
area and degrading the angular resolution. 

The predecessor of the LAT, EGRET [2], which operated between 1991 and 2000, had 
a tracker made of spark chamber modules interleaved with tantalum foils. Its volume was 
80 x 80 x 45 cm 3 and the total on-axis depth was 0.54 radiation lengths, leading to an on-axis 
conversion efficiency of about 30%. Between the tracker and the calorimeter, there was a time 
of flight coincidence system. So they were separated by ~ 80 cm. As a consequence, the field of 
view of EGRET was limited to about ~ 25°. The calorimeter of EGRET was a 77 x 77 x 20 cm 3 
block of Nal, without segmentation. The total on-axis depth was about 8 radiation lengths. Its 
energy resolution was ~ 9% from 100 MeV to 10 GeV. The ACD of EGRET was a monolithic 
dome-shaped scintillator. At high energy, the electromagnetic showers in the calorimeter produce 
backsplash, which can produce some signal in the ACD, which mimics charged cosmic rays. So 
the maximum energy reached by EGRET was limited to about 20 GeV for two reasons: because 
the calorimeter was not segmented (preventing a good energy measurement) and because the 
ACD was not segmented (no ability to reject false vetos from backsplash). 

2.1. The Large Area Telescope 

The LAT was designed to overcome EGRET's limitations. It is composed of 16 identical towers 
arranged in a 4 x 4 array, with each tower made of a tracker module on top of a calorimeter 
module. The overall surface of the LAT is 150 x 150 cm 2 . A tracker module consists of 18 silicon 
(x,y) planes. The first 16 planes are interleaved with a tungsten foil, which thickness is 0.03 
radiation lengths for the first 12 planes and 0.18 radiation lengths for the following 4 planes. 
The total on-axis depth is 1.5 radiation lengths, leading to an on-axis conversion efficiency of 
about 70%. The height of the tracker is 60 cm and there are 5 cm between the tracker and the 
calorimeter. So the overall aspect ratio of the LAT is much smaller than that of EGRET and 
its field of view is about 60°. Each calorimeter module is an hodoscopic array of 8 x 12 CsI(Tl) 
crystals, arranged in 8 layers, each layer being 19.9 mm thick. The total on-axis depth is 8.6 
radiation lengths. Despite of its modest thickness, we will see in section [3] that the calorimeter 
segmentation allows us to measure gamma rays up to 3 TeV with good resolution. The ACD 
is segmented in 89 plastic scintillators tiles in order to minimize its sensitivity to backsplash. 
Thanks to its calorimeter segmentation and its ACD segmentation, the LAT is able to detect 
and measure gamma rays up 3 TeV, much higher than EGRET. 

2.2. The LAT calorimeter 

Each calorimeter module is made of 8 layers with 12 crystals per layer. Each crystal is a 
326 x 26.7 x 19.9 mm 3 block of CsI(Tl). The orientation of the crystals depends on the layer: 



along x for the first layer, along y for the second, and so forth. This geometry allows for full 
imaging of gamma-ray showers. The crystal height is 1.07 radiation lengths and its width is 0.75 
Moliere radii. This thin segmentation gives 2 coordinates of the energy deposition in a crystal. 
The third one, the longitudinal position along the crystal, is given by the asymmetry between 
the signals read out at each end. 

In order to have a large energy range, each crystal end is readout by two PIN photodiodes 
(25 mm 2 and 147 mm 2 ) with two gains per diode, leading to four ranges (2 MeV to 100 MeV, 
2 MeV to 1 GeV, 30 MeV to 7 GeV, 30 MeV to 70 MeV). During normal data taking, the 
lowest non-saturated range is chosen. This very large dynamic range (2 MeV to 70 GeV) allows 
us to measure very low energy deposits (low energy gamma rays, minimum ionizing protons or 
ground muons) while avoiding crystal saturation for gamma rays up to ~ 1 TeV. After 4 years 
of operation, there is only one dead diode out of 6144. 

The in-flight calibration of the crystals is done in several steps. We use a charge-injection 
system to correct for electronics non-linearities. The lowest energy range is calibrated with 
the minimum ionizing cosmic proton peak at 11 MeV. The inter-range calibration is performed 
using non-interacting cosmic protons and heavy nuclei. The latter are also used to calibrate the 
asymmetry between the two ends, with the help of the tracker from which the reconstructed 
direction can be extrapolated to the crystal to give the longitudinal position reference. The 
observed crystal light yield attenuation due to radiation damage is about 1% per year. 




Figure 2. Exploded view of a calorimeter 
Figure 1. Schematic diagram of the Large module showing the hodoscopic 12 x 8 array 
Area Telescope. °f Csl detectors in the carbon fiber structure. 



2.3. Event reconstruction and selection 

The goal of the LAT is to reconstruct and select good photon candidates. We use the tracker 
to find the direction of the event and the conversion point. The energy measurement relies on 
both the tracker and calorimeter. Very low energy gamma rays do not reach the calorimeter 
and their energy is determined from the number of hits in the tracker. For higher energies, we 
correct the energy deposited in the calorimeter for the various losses by using the position of 
the shower maximum or by performing a fit of the longitudinal profile of the shower, as fully 
described in section [3l 

The direction information is used to extrapolate the trajectory of the particle back to the 
ACD in order to reject charge particle background. This information, in addition to the topology 
of the shower in the tracker and in the calorimeter, is used for the selection of the gamma-ray 



candidates. The reconstruction and selection are based on classification trees and have been 
optimized using a detailed GEANT4 simulation. The performance of the instrument has been 
carefully checked in beam tests and on orbit [3] . The on-axis effective area is 0.25 m 2 at 100 MeV, 
0.7 m 2 at 1 GeV and 0.8 m 2 at 100 GeV. 

The rate of cosmic-ray background passing through the LAT can be as high as 10 kHz. Since 
the characteristic signal shaping times of the tracker, calorimeter and ACD are respectivelly 
10, 3.5 and 4 [as, such a rate implies that the probability that particles pass through the LAT 
a few [is before a gamma ray triggers the detector, can be higher than 10%. These out-of-time 
particles can leave sizable signals in all the three subdetectors. These spurious signals can cause 
serious errors in the event analysis and thus must be cleaned by the event reconstruction. In 
the calorimeter this is done by running a clustering algorithm and finding the clusters related 
to the real event. 

3. Energy reconstruction from 3 GeV up to 3 TeV 

The LAT phase space is very large: energies from 20 MeV up to more than 300 GeV and angles 
up to 60 deg with respect to the vertical axis of the detector. Because of the geometry of the 
instrument, the fraction of energy deposited in the calorimeter varies very much throughout this 
phase space. For normally incident photons, the fraction of energy deposited in the calorimeter 
is maximum at 1 GeV and is about 80%. Above 1 GeV, leakage from the calorimeter starts 
to be the main cause of energy loss. Another cause of energy loss is the 4 cm gaps between 
the towers. Here we describe an energy reconstruction method that works for photons above 
~ 3 GeV, when the energy deposited in the tracker is small compared to that in the calorimeter, 
and where the leakage from the calorimeter can be large. 

The longitudinal segmentation of the calorimeter allows a fit of the shower longitudinal profile, 
which is a good way to correct for energy loss when the shower is not fully contained in the 
calorimeter. In order to be able to perform such a fit, we must predict the energy deposited in 
each layer for any longitudinal profile. At normal incidence, the layers are perpendicular to the 
shower axis and the prediction is obvious. But at large incidence angle, the energy contained 
in a slice of the shower at a given age is shared among several layers or part of it escapes the 
calorimeter. As a consequence, for a given shower, we have to model accurately its development 
through the layers of the calorimeter, taking into account the gaps between towers. 

The crystal electronic signals saturate at about 70 GeV. This saturation, which occurs for 
photons above 1 TeV, must be taken into account by the shower fit. When a crystal is saturated, 
we do not know how much energy has been deposited in it, so it cannot be used in the shower 
fit. Since the saturated crystals lie along the core of the shower, the only information we have 
comes from the transverse tails of the shower, which contain less than 10% of the energy. That 
means that, in order to provide a useful energy resolution above 1 TeV, we have to be able 
to model very accurately the development of the shower through the calorimeter at the crystal 
level, and paying close attention to those crystals that lie at more than one Moliere radius from 
the shower axis. 

As a consequence, the shower longitudinal profile fit is based on the modelling of the 
longitudinal and transverse profiles of electromagnetic showers and on the modelling of the 
development of the showers through the LAT calorimeter. 

3.1. Shower longitudinal profile parameterization 

The average longitudinal shower profile can be described by a gamma function [3]: 




(1) 



where t is the longitudinal shower depth in units of radiation length, a the shape parameter 
and (3 the scaling parameter. The position of the maximum of the profile is T = (a — 1)(3. As 
shown in [5], the same function can be used to describe individual profiles. In order to model 
the longitudinal profile, we have to determine how the distributions of a and /3 vary with energy. 

In order to do so, simulations of photons in a very large Csl calorimeter with GEANT4 at 
various energies have been performed. The calorimeter is segmented in 400 layers of 1.85 mm = 
0.1 radiation length. Its total depth is 40 radiation lengths so photons are fully contained. For 
each event, the conversion point is used to define the shower starting point t = 0. In order to 
smooth the longitudinal profile, 0.5 radiation length bins were used to construct the histogram 
of the longitudinal profile, which is fitted with equation JT]). 

Fig [3] shows the distributions of a and j3 and their logarithm for 10 GeV photons. The 
distribution of In a is more gaussian than a, whereas the distribution of j3 is more gaussian than 
ln/3. That is why we decide to use In a and f3. 




Figure 3. Distributions of a and (5 and their Figure 4. Correlation between In a and j3 for 
logarithms for 10 GeV photons. 1,10,100 and 1000 GeV photons. 

Fig [5] shows the correlation of these two variables for various energies. The correlation is 
independent of energy with a slope of tan# c = 0.5. As it is preferable to use uncorrelated 
variables, we use the following variables: 

So = In a cos 6 C + (3 sin 9 C 
S\ = — In a sin 6 C + j3 cos 9 C 

Fig [5] shows the distributions of So and Si for various photon energies. The last step of the 
description of the shower parameters So and Si is the parameterization of the variation of their 
mean and RMS with logE, as shown in Fig [6l The dependence of the mean can be simply 
parameterized with a quadratic polynomial. The dependence of the RMS is quadratic below 
~ 30 GeV and linear above. 

In order to determine the error of the modeling of the shower profile, that is to say the error 
that will be used in the calculation of the x 2 during the fit, we use the residuals of the shower 
profile fit as a function of the shower depth t. The RMS of these residuals are almost flat with a 
broad maximum at 0.5 < t/T < 1.5. As a consequence we assume that the error does not depend 



on the shower depth and we use the same modeling error for all layers. We define €t as the 
RMS of the residuals at shower maximum divided by the energy deposited at shower maximum. 
Again et varies with energy. Its variation can be parameterized by ej<(E) = 0.17e~ 1 °g- E / 1 - 38 ) 
with E in GeV. During the fit, the error used in the x 2 calculation is 5e = €t(E) x max(E'j), 
where Ei is the energy deposited in layer i. 
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Figure 5. Distributions of the shower 
parameters So and Si for 1, 10, 100 and 
1000 GeV photons. 
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Figure 6. Parameterization of the mean and 
RMS of Sq and Si as function of logE. 



3.2. Shower transverse profile parameterization 

The average transverse energy profile can be described by two components, a core and a tail. 
As in [5J, we use the following function: 
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where Rc and Rj< are the medians of the core and the tail components (in units of Moliere 
radius), and p is the relative weight of the core component (0 < p < 1). This transverse profile 
depends on the shower depth t/T. In order to study the dependence of Rc, Rt and p on 
t/T, simulations were performed with GEANT4. In this case, the Csl calorimeter was divided 
longitudinally in slices of 0.1 radiation length and radially in rings with a width of 0.1 Moliere 
radius, which is 3.5 cm for Csl. For each event, we perform the longitudinal profile fit in order 
to determine T, the position of the shower maximum. 

We compute the average fraction of energy contained in a cylinder as a function of the 
cylinder radius for various intervals of the shower depth t/T. Fig. (7|(left) shows some examples 
of these cumulative transverse profiles for 100 GeV photons. The fit of these profiles allows the 
determination of Rc, Rt and p. Fig. [7]^ right) shows how these parameters vary with t/T. The 
variation of the average transverse profile with energy was found to be very weak, so we neglect 
it and use the 100 GeV profile. 



3. 3. Constraining the longitudinal fluctuations during the fit 

The fit of the shower longitudinal profile works well when the shower maximum is well contained 
in the calorimeter. Since the LAT calorimeter is composed of 8 layers, this is the case when 



there is significantly less energy in the last layer than in the seventh layer. For incidence angles 
smaller than 25 deg, the fraction of photons for which the shower maximum is not well contained 
is ~ 25% for photons at 100 GeV. The situation above 1 TeV is even worse since almost no 
photons within 45 deg of the detector axis have a well contained shower maximum. 

In order to be able to reconstruct the energy of these photons, we have to use more information 
than the deposited energies in the layers. The idea is to use the knowledge on the fit parameters 
So and Si (especially the dependence of their mean and RMS with E) by adding to the \ 2 a 
term that constrains the parameters to be close to their expected values. Since Sq and S\ are 
gaussian and uncorrelated, we use Sg + sf, where Si(E) = (Si — /i5 i (S))/o , 1 g.(£?), i.e. the distance 
of Si to its mean value and divided by its RMS. The shower fit \ 2 becomes: 

X 2 (So, S 1 , E)=J2 ^^'J'e)^ 2 + C ( S ° (S) + S " {E) ) 

where e m ^ and e Pi j are respectively the measured and predicted energies in layer i. The parameter 
constraint term is multiplied by a factor c. We varied c and found that c = 1 was optimal. Fig. [8] 
compares the distributions of E^ t /E trvie with c = and c = 1, for 100 and 1000 GeV photons 
at normal incidence. When c = 0, the energy can be overestimated by a factor greater than 2, 
whereas setting c = 1 allows us to remove entirely the high energy tail, even when the shower 
maximum is not well contained. 




Figure 7. Left: average transverse profiles at 
100 GeV. Right: variation of the parameters 
of the transverse profile with the position 
along the shower at 100 GeV. 



Figure 8. Comparison of the distributions 
of Efi t /E tTne with and without the parameter 
constraining term at 100 GeV (left) and 1 TeV 
(right). 



3.4- The development of showers through the LAT calorimeter 

For a given shower longitudinal profile P(E,a,p)-> the predicted energy deposited in layer i is: 

POO 

e P ,i = / fi(t)P(E, a ,p)(t)dt 

J 

where fi(t) is the fraction of energy deposited in layer i by the shower slice between t and t + dt. 
For normal incidence photons, fi(t) only depends on the position of the shower slice t. But when 
the incidence angle is greater than 0, the layers are no longer perpendicular to the axis of the 
shower and so the transverse profile must be taken into account. Since the transverse profile is 
a function of t/T, the energy fraction fi(t) also depends on T. 

The computation of fi(t) depends on the photon trajectory. If the photon converted in the 
tracker, we use the results of the tracking algorithm to get the direction and the conversion point, 
corresponding to the starting point of the shower. We then compute the number of radiation 



lengths crossed by the shower between the conversion point and the extrapolated entry point of 
the photon in the calorimeter. If no track has been found, we use the results of the moments 
analysis of the shower in the calorimeter to get the direction and we assume that the conversion 
point corresponds to the entry point of the photon in the calorimeter. 

Using this trajectory information, the computation of fi(t) is done as follows (illustrated in 
Fig. ED: 

• divide the trajectory into ~ 1.85 mm steps, starting 100 cm ahead of the entry point of the 
photon in the calorimeter; 

• at each step along the trajectory, at a distance x in mm along the trajectory from the 
starting point of the shower, corresponding to the shower depth t: 

— define a disk centered on the step point, perpendicular to the trajectory, whose radius 
is 3 Moliere radii = 95 mm (the red line in Fig. EJ ; 

— define a grid of test points on the disk, with a step of 0.05 Moliere radius; 

— the test points are used to compute the fraction of energy deposited in the layers at 
the shower depth t, using the transverse profile parameterization depending on t/T 
(represented in Fig. El by the red area); 

— compute di, the number of radiation lengths corresponding to the longitudinal step, 
taking into account the active (Csl cystals) and passive (gaps between calorimeter 
modules) material seen by the shower, and add it to t. 

While developing through the calorimeter, the shower becomes wider. At some point, two 
different parts of the shower at the same shower depth t evolve independently. In other words, 
if one part of the shower escapes the calorimeter, the rest of the shower evolves as if the whole 
shower was inside the calorimeter. This is taken into account when computing dt. At this stage 
of the computation, we have the energy fractions as a function of x as well as the shower depth 
t as a function of x. We then deduce the energy fractions as a function of t: fi(t). 

Examples of fi(t), the fractions of energy deposited in the layers, are shown in Fig. \10\ for a 
shower maximum T equal to 8 radiation lengths. The fraction of energy deposited in layer i at 
normal incidence is very simple, being 1 when the shower is inside layer i, and otherwise. For 
an incidence angle of 45 degrees, the situation is more complex, because of the sharing of the 
energy between layers, which increases with t as the shower becomes wider. 




Figure 9. One step of the computation 
of the development of an individual 
shower inside a calorimeter module. 
The arrow corresponds to the photon 
trajectory. The red disk represents a 
shower slice at a given shower depth for 
which we compute the fraction of energy 
deposited in each layer. 




Figure 10. Predicted fraction of energy deposited 
in the layers of the calorimeter for a photon at 0° 
(left) and 45° (right) cLS Si function of the shower 
depth. Each curve corresponds to one layer: from 
the first layer on the left to the last layer on the 
right. 



3.5. Crystal saturation handling 

Due to the readout limitation, the signal of a crystal saturates when the deposited energy is 
greater than ~ 70 GeV. It occurs for photons with E > 1 TeV. In order to take this saturation 
effect into account in the shower fit, we consider the saturated crystals separately, and especially 
take out their energy from the layer energies. The computation of the shower development 
through the calorimeter (presented in section I3.4p is used to determine the fraction of energy 
deposited in the saturated crystals as a function of t. The standard part of the \ 2 becomes: 

X 2 = ( em >* ~ e P^ 2 / 6e2 + Y (max(0, e m ,j - e P)i )) 2 /5e 2 

layers sat. crystals 

where the layer energies are summed over the non saturated crystals and where we introduced 
the term max(0, e m ^ — e Pi j) in order that the saturated crystal contribution vanishes as soon as 
the predicted energy is greater than the measured energy. 

3.6. Results 

The performance of the shower profile fit is summarized in Fig [11] The bias is less than 3% 
and is easily corrected for. The energy resolution varies with energy and incidence angle. At 
normal incidence, it degrades when the energy increases. This is because the shower containment 
decreases with energy, as shown in Fig 1 1 1 If right ^ . At 45°, the resolution improves slightly from 
3 GeV to 30 GeV because the fraction of energy deposited in the tracker decreases. Above 
30 GeV, the resolution starts to degrade because the shower containement decreases. It must be 
noted that at 3 TeV, the shower containement is only ~ 10% and the mean number of saturated 
crystals is ~ 6. But still we have a decent energy resolution, especially for off-axis photons. 
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Figure 11. Performance of the shower profile fit: bias (left), resolution (center), shower 
containment (right). 



4. Some Fermi results at high energy 

It is impossible to summarize all the results obtained by Fermi after 4 years of operation. We 
only focus on results at high energy for which the measurement of the energy is important. 

4-1. In-flight absolute energy scale 

There is no gamma-ray source with a sharp and well known spectral feature for absolute energy 
scale calibration. But we can use the cutoff of the electron spectrum due to the Earth's magnetic 
field to check our absolute energy scale [6J. As explained in section l4To] Fermi can detect electrons 



and positrons as efficiently as gamma rays. The geomagnetic cutoff rigidity depends on the 
geographic longitude and latitude. The altitude of the Fermi orbit is 565 km and its inclination 
is 26.5°. So the Fermi spacecraft spans a range of geomagnetic cutoff that goes from 6 to 14 GeV. 
This range is divided into 6 bins. In each geomagnetic cutoff bin, we first determine the expected 
primary electron+positron spectrum. This is done by using a tracer code that computes the 
trajectories in the reverse direction (from the spacecracft to space) in order to reject secondary 
particles that are produced in the atmosphere. The measured primary spectrum is obtained 
after subtraction of the secondary fraction which is derived using the azimuthal distribution. 
The ratio of the cutoffs of the expected and measured spectra is shown in Fig. [TJ] From this 
comparison, we find that, in the energy range [6,14 GeV], the energy scale is known within 5%. 

4-2. Gamma-ray sources spectral features 

By detecting more than 1800 sources, Fermi is able to peform detailed population studies. 
Concerning active galactic nuclei, the most numerous sources for Fermi, spectral breaks have 
been observed to be common in Flat Spectrum Radio Quasars (the break energy ranging from 
1 to 10 GeV in the source frame for the brightest sources), and present also in some bright 
low-synchroton peaked BL Lacs [7). These breaks represent challenges for theoretical models 
aiming at describing the blazar phenomenon. Another important class of objects are pulsars: 
Fermi has currently detected more than 100 pulsars, whereas EGRET only detected 6. Pulsars 
have a power-law spectrum with a cutoff, which position and shape depends on the location of 
the gamma-ray emission region in the pulsar magnetosphere. The cutoffs of all pulsars measured 
by Fermi [8] are consistent with exponential cutoffs at a couple of GeV. This result rules out 
the low-altitude emission models which predict a sharper cutoffs. 

4-3. Gamma-ray bursts and Lorentz invariance violation 

Quantum gravity may cause the violation of Lorentz invariance and the speed of light could then 
depend on energy. It would imply that photons with different energies emitted simultaneously 
by a very distant source would not arrive at the same time. Fermi is able to constrain the 
violation of Lorentz invariance by looking at the highest energy events of gamma-ray bursts. 
GRB090510 [9] was a short gamma-ray burst with a redshift z = 0.9. The highest energy event 
has 31 GeV and is detected 0.86 ms after the start of the emission below 1 MeV (measured by 
the Gamma-ray Burst Monitor, the second instrument on-board Fermi). In the framework of 
linear Lorentz invariance violation, the delay between two photons is proportional to the energy 
difference divided by Mqq, the quantum gravity mass scale. Assuming that the high energy 
emission did not start before the burst, this 31 GeV event sets the limit: MQQ,/Mpi anc k > 1.2. 

4-4- Earth limb spectrum 

Gamma rays from the Earth limb are the result of the interaction of cosmic rays in the 
atmosphere. Above 10 GeV, the emission that is detected by Fermi is only due to grazing 
incidence cosmic rays in the upper atmosphere. The altitude range of the emission is rather 
thin: cosmic-rays passing too high do not interact enough and the gamma-rays produced by 
cosmic-rays that pass too low are attenuated by interacting in the atmosphere. The measured 
Earth limb spectrum [10] is a power law with spectral index = 2.79 ± 0.06, in agreement with 
the spectral index of cosmic-ray proton = 2.75. The spectrum is very smooth up to 1 TeV. The 
absence of any significant sharp feature is an additional proof that our photon selection is well 
under control. 

4-5. Extragalactic diffuse emission 

The largest source of Fermi photons is the galactic diffuse emission, due to the interaction 
of cosmic-rays with the matter in the Galaxy. But there is a much fainter diffuse emission 



that has an extragalactic origin and that represents the history of the non-thermal universe. 
Thanks to a careful background subtraction, we measured the extragalactic diffuse emission 
up to 100 GeV and expanding this measurement up to 1 TeV is under way. Currently, 
unresolved sources like blazars, radio galaxies and star-forming regions account for 50% to 80% 
of this emission. 

4-6. Cosmic-ray electron and positron spectrum 

Since Fermi's primary goal is to detect gamma-rays, the main trigger configuration rejects a 
large fraction of charged particles by requiring that there is no activity in the ACD above the 
tracker tower that triggered. But, in order to be fully efficient at high energy, all events that 
deposit more than 20 GeV in the calorimeter are sent to the ground, which is the case for high 
energy electrons and positrons. Since they produce electromagnetic showers very similar to 
photon induced showers, Fermi is able to reconstruct and select electrons and positrons a well 
as gamma-rays. The fraction of residual background is estimated to vary smoothly from 4% at 
20 GeV to 20% at 1 TeV. Since measuring such high energies with a 8.6 radiation length thick 
calorimeter is not an easy task, we checked that our results did not change when selecting off- 
axis events that pass through more than 12 radiation lengths in the calorimeter. The measured 
spectrum [12] rules out the excess around 700 GeV seen by the ATIC experiment [13] and can 
be explained by a conventional model with an adjusted injection spectrum or with an additional 
source, which could be a nearby pulsar. 
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Figure 12. Distribution of cutoff energy 
ratios as a function of Mcllwain L (with 
corresponding cutoff energies indicated). The 
statistical errors are represented by the black 
error bars, and the systematic errors are in 
gray. From [6j 
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Figure 13. Cosmic-ray electron spectrum 
as measured by Fermi LAT for 1 yr of 
observations along with other recent high- 
energy results. Systematic errors are shown 
by the gray band. From [12] . 



4-7. Cosmic-ray positron fraction 

Without any magnet on-board, Fermi is not able to discriminate between electrons and positrons. 
But because of the geomagnetic field, there is a region of the Earth's horizon from which only 
electrons can come and another region from which only positrons can come, as shown in Fig. O 
So Fermi is able to measure separately the spectra of electrons and positrons between 20 GeV and 
200 GeV (above 200 GeV these regions become too small). From these spectra we can deduce 
the positron fraction, which is shown in Fig. [151 Our measurement [14j . which was the first 
above 100 GeV, confirms the increase of the fraction with energy as reported by PAMELA |15j . 
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Figure 14. Examples of calculated electron 
(red) and positron (blue) trajectories arriving 
at the detector, for 28 GeV particles arriving 
within the Equatorial plane (viewed from 
the North pole). Forbidden trajectories are 
solid and allowed trajectories are dashed. 
From [T4] . 
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Figure 15. Positron fraction measured by 
the Fermi LAT and by other experiments. 
The Fermi statistical uncertainty is shown 
with error bars and the total (statistical plus 
systematic uncertainty) is shown as a shaded 
band. From |14j . 



5. Conclusion 

After 4 years in orbit, Fermi has proved to be an excellent gamma-ray telescope. Thanks 
to its large effective area and field of view, to its good angular and energy resolution, Fermi 
has already unravelled the gamma-ray sky with much more precision than its predecessor 
EGRET. By detecting more than 1800 sources, it allows population studies of pulsars, supernova 
remnants, active galactic nuclei and gamma-ray bursts. The diffuse emission of the Galaxy and 
the extragalctic diffuse emission are also known much better. Fermi has proved to be a very 
good detector of electrons and positrons as well, and, by using the geomagnetic field, it is even 
able to discriminate between them. Part of the success of Fermi is due to its large energy range. 
Thanks to the fine segmentation of its calorimeter, we are able to measure and select gamma 
rays, electrons and positrons up to 3 TeV, despite the modest thickness of the calorimeter. 
Taking benefit from our experience after 4 years in orbit, we are currently optimizing the event 
reconstruction and selection. With a better sensitivity, Fermi will certainly continue to improve 
our understanding of the high energy sky. 
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